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System studies show that a N3-X hybrid wing-body aircraft with a turboelectric 
distributed propulsion system using a mail-slot inlet/nozzle nacelle can meet the 
environmental and performance goals for N+3 generation transports (three generations 
beyond the current air transport technology level) set by NASA’s Subsonic Fixed Wing 
Project. In this study, a Navier-Stokes flow simulation of N3-X on hybrid unstructured 
meshes was conducted, including the mail-slot propulsor. The geometry of the mail-slot 
propulsor was generated by a CAD (Computer-Aided Design)-free shape parameterization. 
A novel body force model generation approach was suggested for a more realistic and 
efficient simulation of the flow turning, pressure rise and loss effects of the fan blades and 
the inlet-fan interactions. Flow simulation results of the N3-X demonstrates the validity of 
the present approach. An optimal Shape design of the mail-slot nacelle surface was 
conducted to reduce strength of shock waves and flow separations on the cowl surface. 
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Nomenclature 

sectional area normal to the axial direction 

number of blades in a row 

boundary layer ingestion 

Computer Aided Design 

blade chord length 

body force 

adjustment function in spanwise direction for normal body force coefficient (see Fig. 3) 
adjustment function in local mass flow rate for normal body force coefficient (see Fig. 5) 
height 

hybrid wing-body aircraft 
blade-to-blade gap- staggered spacing 
body force coefficient 
length 

Mach number 

mass flow rate 

meridional coordinate 

pressure 

gas constant 

radius 

entropy 

temperature 

velocity 

cylindrical coordinates 
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a = blade camber angle 

y = ratio of specific heats 

p = density 

a = solidity 


Subscripts 
c - 

h - 

LE 

n : 

p 

rel - 

TE 

t 

tot - 


corrected 
hub of blade 
leading edge 

direction normal to local flow in a blade-to-blade section 

direction parallel to local flow in a blade-to-blade section 

relative to the blade 

trailing edge 

tip of blade 

total 


I. Introduction 

T he Subsonic Fixed Wing Project of the NASA Fundamental Aeronautics Program has focused on research of 
the next generations of civil aircraft. For the near term (N+l), midterm (N+2), and far term (N+3), specific 
goals have been defined for reducing noise, emissions, and fuel burn from the current state-of-the-art aviation 
technology. 1 To achieve our aggressive goals, especially for the N+2 and N+3 time frames, we need to explore 
innovative design concepts in airframes and propulsion. 

The hybrid wing-body (HWB) aircraft, although not new, is an alternative design concept to the conventional 
tube-and-wing aircraft. 2 The fuselage and wings of the HWB are integrated into a flying wing, which has better 
aerodynamic efficiency than tube-and-wing configurations. Also, the wide airframe body of the HWB configuration 
is beneficial for shielding downward-propagating noise coming from the engines that are installed above the aircraft. 

In comparison to conventional pylon-mounted engines, embedded engines with boundary-layer-ingestion- (BLI) 
offset inlets reduce ram drag, wetted area, structural weight, and noise. 3,4 Disadvantages of the embedded engines 
are higher flow distortion and reduced pressure recovery at the engine faces because of the BLI-offset inlet. 

The N3-X was proposed as a candidate to meet the N+3 goals. 5 The N3-X is a 300-passenger HWB aircraft 
employing turboelectric distributed propulsion (TeDP), which utilizes superconducting electric generators, motors, 
and transmission lines (see Fig. 1(a)). The TeDP system allows power generation and thrust generation to be 
separated and enables a small number of turboshaft engines to drive any number of propulsors, each of which is 
composed of an electric motor and a fan. The N3-X was designed to have two large turboshaft engines at each wing 
tip to feed the engines with uniform free-stream flows. In addition, BLI propulsors are distributed in the two- 
dimensional mail- slot inlet/nozzle nacelle on the upper surface of the HWB near the trailing edge. A section view of 
an electric-motor-driven propulsor is depicted in Fig. 1(b). According to the system study in Ref. 5, the N3-X is 
expected to reduce mission fuel burn by more than 70 percent relative to a reference aircraft, a Boeing B777-200LR 
flying 7500 nm with a cruise of Mach 0.84 and a 118,100-lb payload. 



(a) Overall configuration. (b) Section view of nacelle including fan and electric motor. 

Figure 1. N3-X configuration with turboelectric distributed propulsion. 5 
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The conventional uniform-back-pressure boundary condition for fan faces has been used widely to simulate flow 
for propulsion-airframe integration problems. 6 ' 8 However, for integration problems with highly distorted onset flows 
at the fan faces, the assumption of uniform static pressure at the fan faces may not be valid. The low momentum at 
the low-recovery region at a fan face causes a larger local incidence angle for the fan rotor blade. This causes the fan 
to have a stronger suction effect, which, in turn, mitigates the flow distortion. Consequently, use of the uniform 
back-pressure boundary condition has produced conservative results for the flow distortion at fan faces, although it 
is still valid for the qualitative evaluation of inlet fan interaction problems when the diffuser flow separation does 
not reach to the fan face. 9 

A direct coupling 10 of the inlet and full-annulus fan blades in the computational domain would give a more 
realistic simulation of the inlet-fan interaction, but the computational cost would be prohibitively large, especially 
for shape-design applications, in which more than tens of flow simulations are usually required. The body force 
approach 9,11 ' 14 has been drawing attention as an alternative to simulating the coupling with full-annulus fan blades. 
This approach uses body force terms to model flow turning and loss due to rotor/stator blade rows. The body force 
terms are added as source terms in the flow equations for grid cells swept by blade rows. Body force coefficients or 
parameters need to be fitted to, or interpolated from, single-passage Navier-Stokes flow simulation results or 
experimental test data. The body force approach allows relatively accurate flow simulation of BLI inlet-fan 
interaction problems that consider blade force effects without actual full-annulus simulation of the rotor/stator blade 
rows. 

In the present study, we conducted high-fidelity flow simulations of the N3-X configuration to understand the 
effects of the propulsion-airframe integration on the flow field and performance of the vehicle and the distributed 
propulsion system. A CAD-free shape parameterization was used to efficiently generate and modify the complex 
mail-slot nacelle geometry. An efficient body force approach was suggested to incorporate the fan blade effects of 
turning and flow loss. Then the nacelle cowl shape was designed to improve aerodynamic performance of the 
nacelle. 

The remainder of this paper is organized as follows. Section II presents numerical approaches for simulating 
flow, including flow solvers, mesh generation method, and the body force method. Section III describes formulation 
of the present body force approach. Section IV presents validation results of the present body force method. Section 
V explains shape parameterization and installation approaches for the mail-slot nacelle. Section VI presents flow 
simulation and shape design results for the N3-X configuration. Finally, Section VII provides the summary and 
conclusions. 


II. Methodologies for Flow Simulation 
A. Unstructured-Grid Flow Simulation 

GO -flow, 15 a finite- volume unstructured-grid Navier-Stokes solver, was used for all the flow simulations in this 
study except turbomachinery single passage simulations. The compressible Reynolds-averaged Navier-Stokes 
equations were discretized by the cell-vertex finite volume method. Control volumes were nonoverlapping dual cells 
constructed around each node. Each edge connecting two nodes was associated with an area vector of the control 
surface, at which flow fluxes were computed. For second-order accuracy, a linear reconstruction of the primitive 
gas-dynamic variables inside the control volume was used in conjunction with a limiter. The inviscid flux was 
computed using approximate Riemann solvers, and Menter’s two-equation Shear Stress Transport (SST) model was 
used to incorporate turbulence effects. 16 For the time integration, the Lower-Upper Symmetric Gauss Seidel (LU- 
SGS) implicit method was adopted. Parallel processing was accomplished by domain decomposition of the 
computational mesh and the Message-Passing Interface (MPI). 

An unstructured computational mesh was generated with the Mixed-Element Grid Generator in 3 Dimensions 
(MEGG3D) code. 17,18 A triangular surface mesh was generated directly on a background unstructured mesh or from 
stereolithography (STL) data, which can be exported from a CAD model. Then, advancing front/layer methods were 
used to make a hybrid volume mesh. Layers of prismatic cells were generated near viscous walls, tetrahedral cells 
were generated in the remaining computational domain, and pyramid cells were generated between the other cells 
when necessary. 


B. Single-passage Tubomachinery Flow Simulation 

For single-passage Navier-Stokes flow simulations of the R4 fan, the SWIFT code 19,20 was used. The SWIFT 
code is a three-dimensional, multiblock, structured-grid Navier-Stokes analysis code for turbomachinery blade rows. 
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The code solves the Navier-Stokes equations on body-fitted grids using an explicit finite -difference scheme. The 
AUSM+ scheme 21 was used for spatial discretization of the inviscid flux. For the thin-layer approximation, viscous 
terms were included in the blade-to-blade and hub-to-tip directions and were neglected in the streamwise direction. 
For turbulence effects, Wilcox’s 2006 k-co model with a stress limiter was selected. The discretized equations were 
solved with a multistage Runge-Kutta scheme with local time stepping and residual smoothing being used to 
accelerate convergence. 


III. Body Force Model 

The body force model used for this study was based on the model developed by Gong, 12 which makes the 
assumption of an infinite number of blades and axisymmetric flow in each infinitesimal blade passage. The cascade 
blade forces in the relative frame of reference are modeled as normal and tangential forces to the local flow, which 
are similar to the lift and drag forces of an airfoil. 


A. Normal Force Model 

In the present study, a modified formulation in Ref. 22 was used to model the normal force: 


K 


( Aa 


F n = — VV +-sin — 
n h n p c { 2 


V 2 


( 1 ) 


where K n is the normal force coefficient. Equation (1) is similar to Gong’s formulation except for the second term of 
the right-hand side. In Eq. (1), h is the blade-to-blade gap-staggered spacing given by 


27Ir^/a 


cos a 


B 


( 2 ) 


where r is the radius, a is the solidity, a is the local blade camber angle, and B is the number of blades. In Eq. (1), 
V n and V p are velocity components normal and parallel to the local cascade flow and are defined as follows: 

V n = Vq cos a - V x sin a (3) 

V p = V x cos a - Vq sin a (4) 

where V x and V e are velocity components in the axial and circumferential directions. In the second term on the 

right-hand side of Eq. (1), c is the chord length and Aa is the camber angle difference between the trailing edge TE 

and leading edge LE\ 


Aa — clte ~ ^le 


( 5 ) 



Figure 2. Illustration of body force components at blade-to-blade section 
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Axial and circumferential components of the normal force are calculated by 


CD IS 

* 1 ^ 

II 

(6) 

II 

(7) 


where V rei is magnitude of the blade relative velocity. 


Defoe 23 suggested an empirical model for determining the normal force coefficient K n for a particular fan rotor: 


K n = (4.2- 3.3a j 4.172 


r-n 


-3.118 


r-r,. 


+ 2.145 


( 8 ) 


where r h and r t are the hub and tip radius of the fan blade, respectively. The first expression in parentheses in 
Eq. (8), an empirically obtained term suggested in Ref. 12, is multiplied by the bracketed expression to adjust the 
magnitude of K n along the spanwise direction. For this study, we adopted a formulation for K n that is very similar to 
Eq. (8) but allows more flexibility: 


K n = (4.2 - 3.3a)/(r) (9) 

where /(r) is a set of line segments connecting control points distributed along the blade span as shown in Fig. 3. The 
control points of /(r) were adjusted to match available computational fluid dynamics (CFD) or experimental data as 
will be shown in Section IV A. 



Figure 3. Definition of f(r) for the body force coefficient K n . 


B. Parallel Force Model 

In Gong’s model, F p , the component of the body force for the flow loss, is obtained by 



( 10 ) 


with K p = 0.05. Marble 24 suggested a loss force model relating the entropy production and the parallel body force as 
follows: 


_ t Yjl 
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where T is the temperature, s is the entropy, and m is the coordinate along the meridional streamline and V m is 
its velocity component. In this study, Marble’s parallel body force model was employed instead of Gong’s model 
because the direct relation between the parallel body force and entropy production in Eq.(ll) allows an easier and 
physically more meaningful modeling process. 

In Refs. 10 and 11, span wise distributions of the entropy production across the blade region were used to 
determine the flow loss. Here, we made an assumption that the entropy production is constant along the blade span 
and chord for a simpler modeling procedure. Thus, we could calculate the averaged entropy change, either from test 
data or CFD results for the total pressure P tot and total temperature T tot ratio of a blade row, by using the canonical 
equation of state: 


A/? 

R 


y - 1 


-In 


T 

v 1 tot\ y 


-In 


( 12 ) 


where R is the gas constant, y is the specific heat ratio, and stations 1 and 2 refer to the entrance and exit of a fan 
blade row, respectively. As an example, Fig. 4 shows entropy production values across a fan rotor blade row 
calculated by Eq. (12) with the total pressure and total temperature ratios from test data. 25 In a fixed-rotor- speed line, 
the entropy production decreases until the minimum loss point and then increases as the MFR increases. 



Figure 4. Mass flow rate versus entropy production of the R4 fan rotor 
blade row for different rotor speeds. 


The entropy gradient in Eq. (11) is approximated by the entropy production of a blade row and the section chord 
length in the meridional plane at each radial location: 


ds ~ As 

dm x TE -x LE (] 3 ) 


The three components of F p are calculated as follows: 


F =F 

P,X ± p 


Vx_ 


F p,q = F P 


Yl 
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and 


F =F 

P,r ± p 


Yl. 


(14) 


The sum of entropy production of a blade row changes as the MFR is varied at a constant rotor speed. An 
entropy production curve is a function of the corrected MFR. To make actual interpolations along the distributions, 
we need an interpolation parameter. 

For structured-grid approaches having axisymmetric grid sections, each grid section is treated as an independent 
axisymmetric flow section, and the sectional MFR is used as the interpolation parameter for cases with a distorted 
inflow. 11 However, validity of the parametric use of the sectional MFR for the entropy production is questionable. 
The use of the sectional average MFR would cause a smaller entropy production than the actual sum of local entropy 
production at each point of a section due to the convex nature of the entropy production distribution shown in Fig. 4 
especially for a distorted inflow condition. Therefore, in this study, we defined a local MFR parameter as 


"hoaa=pvAx) 


( 15 ) 
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where p and V x are the local density and local axial velocity, respectively, and A(x) is the flow passage area normal 
to the axial direction. ril local is also corrected in the same way as the global mass flow rate. 

An alternative interpolation parameter available to relate the body force model and local flow conditions is the 
local flow coefficient 14 defined by a local axial speed divided by a local rotor speed. The present local MFR 
parameter allows a direct usage of entropy production curves such as shown in Fig. 4, and thus a body force model 
can be generated in a straightforward way. 

The normal body force coefficient K n is also augmented by multiplying an adjusting function g( m local ) shown in 
Fig-5. 


K n = (4-2 - 33d)f(r)g(m local ) 


(16) 



Figure 5. Illustration of an adjusting function g( m local ) for a normal body force model. 

C. Summary of the Procedure 

A summary of the procedure for generating the body force model follows: 

1) Obtain an entropy production curve for a given fan geometry and rotor speed from test data or single - 
passage Navier-Stokes analyses as shown in Fig. 4. 

2) Set an initial entropy production model with asymptotic lines for the target entropy production curve as 
shown in Fig. 6. Set entropy production ds = max(dsl,ds2) at each m local - 

3) Iteratively run a three-dimensional Euler solver with the body force model and adjust slopes and intercepts 
of the entropy model lines to match the resulting entropy changes to the target entropy production 
distribution. 

4) Iteratively run the Euler solver with the body force model to adjust the control parameters of function 
g(m local ) i n Eq. (16) to match the fan pressure ratio for the speed line. Adjust each control point of f{r) 

separately if spanwise swirl and total pressure distributions are available and need to be matched to each 
other. Adjustments of fan pressure ratio parameters may require additional iterations on the entropy 
production model. In that case, go to Step 3). 

5) If the body force model will be used in Navier-Stokes simulations, run a few Navier-Stokes simulations 
with the body force model to quantify the effects of the hub and shroud boundary layers on the resulting 
entropy production. Pull down the entropy curve functions to compensate for the redundant entropy 
increment in the Navier-Stokes analysis. 



Figure 6. Illustration of initial model for entropy production variation w.r.t. local mass flow rate. 
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IV. Validation of the Body Force Method 

A. Uniform Inflow case 

The R4 fan stage, 25 which was designed by General Electric Aircraft Engines and tested by the NASA Glenn 
Research Center, was selected to validate the present body force model. The R4 fan is a one-fifth scale-model 
representation of the bypass stage of a current-generation high-bypass turbofan engine. The fan has a 22-in. diameter 
and 22 wide-chord rotor blades. A Sectional side view of the R4 fan and nacelle is shown in Fig. 7 (a). The stage 
design-point pressure ratio was 1.47 at a corrected speed of 12,657 rpm.. 

For single-passage Navier-Stokes flow simulations of the R4 fan, the SWIFT code 19,20 was used. The SWIFT 
code is a three-dimensional, multiblock, structured-grid Navier-Stokes analysis code for turbomachinery blade rows 
and described earlier in Section II B. 

An unstructured grid with tetrahedral cells was generated around the nacelle cowl, shroud, and hub of the R4 fan 
test geometry, and GO -flow was run in an inviscid flow simulation mode with the body-force source terms added. 
For the fan downstream exit boundary condition, static back pressure was imposed with the radial equilibrium 
condition. Figure 7(b) shows the computational mesh for the clean-inflow R4 case, and Fig. 7 shows the total 
pressure contours normalized by the free-stream total pressure. The total pressure increases along the streamwise 
direction with an almost uniform slope. 

Figure 8(a) compares the fan rotor pressure ratio versus the corrected MFR at 85 percent of the design rotor 
speed. The test data in Ref. [25] agree very well with the results of the Swift code and the inviscid simulation with 
the body force model. Figure 8(b) shows that the entropy production calculated by the Euler + body force model 
matches the test data well. The averaged entropy production of the R4 fan rotor decreases as the MFR increases until 
the peak efficiency point — approximately 42 kg/s. For MFRs higher than the peak-efficiency MFR, entropy 
production is proportional to the MFR because the incidence angle of the blade section deviates from the minimum 
loss angle. 




Figure 6. Geometry and surface mesh of the R4 fan and nacelle. 
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Figure 7. Euler simulation of R4 fan + nacelle with body force model: 
Normalized total pressure contours at a center section. 



(a) Mass flow rate versus fan pressure ratio for 85 -percent speed. 



(b) Mass flow rate versus entropy production. 


Figure 8. Comparison of body force model results with test data 25 and the Swift code at 85-percent speed. 
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In Fig. 9, spanwise distributions of flow characteristics circumferentially averaged at the rotor exit plane are 
compared for the Swift code and the present Euler + body force model. Actually, the control points of /(r) in 
Section III A are adjusted to closely follow the trend of the swirl and the total pressure distributions for the flow 
conditions near the peak efficiency (m c = 43kg/s). Therefore, the two results are in good agreement except near the 

hub and shroud walls because inviscid simulations were used with the body force approach for a fast turnaround 
time. 

As the corrected MFR decreased to 38 kg/s, as shown in Fig. 10, deviations between the Swift code and 
Euler+body force results became larger, but the trend was still well captured by the body force model approach. If 
the control points of /(r) were also varied w.r.t the local mass flow rate, the body force model approach results 
would be improved even further. The adjustment function g( m local ) was not applied in this case. 



(a) Swirl-angle distribution. (b) Total pressure ratio distribution. 



(c) Total temperature ratio distribution. 

Figure 9. Comparison of rotor exit flow characteristics 
along the span near the choking condition ( rh c =43 

kg/s). 
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(c) Total temperature ratio distribution. 

Figure 10. Comparison of rotor exit flow characteristics 
along the span near the surge condition ( rh c =38 kg/s). 


B. Non-uniform inflow case: Inlet A + R4 fan 

In Ref. [10], an unsteady unstructured grid approach was applied to simulation of a combined inlet-fan system 
including the whole annulus of blade rows. The mesh size was about 59.1 million. The rotor speed was 100% design 
speed of the R4 fan. The S-duct inlet A was scaled up and connected to the R4 fan stage as shown in Fig.l 1 

Their results qualitatively shows that the fan performance was degraded as the distorted flow through the S-duct 
is fed to the fan. This case provides an excellent testbed to the validation of the body force method for distorted 
inflow conditions. 

As a first step for the validation of the present body force method, a body force model was generated for the 
100% design rotor speed of the R4 fan stage. Figure 12 shows the fan pressure ratio and entropy production 
variations are well matched both for fan rotor and stator stages. With the generated body force model, the inlet A + 
R4 fan stage simulations were conducted to see the validity of the model for the distorted inflows. Mesh size for the 
present simulation was about 30 million. Figure 13 compares the fan pressure ratio versus mass flow rate for clean 
and distorted inflow conditions. The present body force method shows very good agreement with the full annulus 
CFD results in Ref. [10]. In Fig. 14, fan adiabatic efficiency was compared between the distorted inflow results by the 
body force approach and the clean inflow results by the test 25 . The distortion by the inlet A caused about 4% point 
reduction in the fan efficiency. 
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Figure 15 compares mach contours at center plane. The results by the body force method gives qualitatively 
good agreements with the full annulus CFD results in Ref. [10]. Figure 16 compares Mach contours at the rotor-stator 
interface location (the vertical blue line in Fig. 11) by the full annulus CFD 10 and the present body force method. 
The full annulus CFD captures rotor blade wakes because it performs unsteady computation. Other than that, overall 
trend of the distorted Mach contours are qualitatively in good agreement. 



Figure 11. Side view of the combined system of the R4 fan stage and scaled inlet- A. 10 



(a) Mass flow rate versus fan pressure ratio. (b) Mass flow rate versus Entropy production. 


Figure 12. Comparison of body force model results with test data at 100-percent speed with clean inflow . 
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Figure 13. Comparison of body force model results with test data and full annulus CFD results 
at 100-percent speed for clean and distorted inflow conditions. 



Figure 14. Comparison of stage efficiency for clean and distorted inflow conditions 

at 100-percent speed. 


13 


This material is declared a work of the U.S. Government and is not subject to copyright protection in the United States. 


Downloaded by NASA GLENN RESEARCH CENTER on July 22, 2013 I http://arc.aiaa.org I DOI: 10.2514/6.2013-2413 



(d) Near stall: present body force method. 


Figure 15. Comparison of Mach number contours 
by full annulus CFD 10 and the present body force method at center section 
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Figure 16. Comparison of Mach number contours at rotor-stator interface position 


V. Geometry Generation of the Mail-Slot Nacelle for N3X 

A Fortran code was written that uses shape parameters to build mail -slot nacelles. The mail-slot nacelle shape 
can be divided into outer cowl and flow passages, each of which is composed of an S-duct inlet, a nozzle, a 
cylindrical duct connected to the inlet and nozzle, and a motor fairing center body. The S-duct, which is a 
transitional duct, has a sectional shape that transitions from a rectangle to a circle. The nozzle, which is also a 
transitional duct, has a sectional shape that transitions from a circle to a rectangle. In both the inlet and nozzle, the 
flow passages are separated by splitters. 

The mail-slot propulsor is built up by 206 surface patches. Each surface patch has a structured plot3D data 
format. The neighboring patches are precisely joined without any gaps so that they are “water tight” and ready for 
unstructured surface grid generation. Figure 17 shows some parameters for the mail-slot propulsor geometry. Other 
parameters such as sectional areas also are defined. Any new parameters can be readily added if necessary. The 
computer code can be used for parametric design and optimal shape -design studies for distributed propulsion 
systems with a mail-slot nacelle. In the present study, the mail-slot nacelle has 16 fans and is swept forward by 25 
degrees. 

Once the water-tight nacelle surface patches are generated, the nacelle surface meshes are combined with an 
airframe surface mesh for a clean HWB configuration. After this installation, a computational unstructured surface 
mesh is generated on the combined configuration. The HWB airframe can be parameterized independently from the 
propulsor then combined with the propulsor to build the full N3-X configuration. 


^ nacelle 



Figure 17. Side section view of schematic representation of the propulsor parameterization. 
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Figure 18 shows the HWB configuration with the distributed propulsion nacelle installed on the upper surface. 

In the present study, the installation was made so that the leading edge of the mail-slot nacelle is located at the 85% 
chord position at the airframe center section. The streamwise location of the nacelle also can be modified by 
changing the parameter xinlet_LE in Fig. 17. The nacelle cowl follows the span wise variation of the wing upper 
surface. Also, the cowl lips are twisted down in the outboard direction according to local slopes at the inlet leading- 
edge on the wing surface. Figure 1 8 (d) shows the computational mesh generated for the combined configuration. 
About 31 million computational mesh points are in the half model of the N3-X configuration, and the first nodes off 
the viscous walls are clustered to the wall so that the y+ values at the first nodes are less than 2. At this point, the 
flow simulation could be conducted with the proper body force model for the fan blades. 



(a) Mail slot nacelle installed on the JN3-X. (b) Inside the flow passages. 



Figure 18. N3-X configuration with installed mail-slot propulsor 


VI. Flow Simulation and Shape Design Results of N3-X 
A. Body Force Model Generation 

The overall size and performance specifications of the propulsor fan for N3-X were interpolated from data 
presented in Table 5 of Ref. 5: fan diameter = 1.02235 m (40.25 in.), rh c = 150.8 kg/s, pressure ratio = 1.325, and 

rotor speed = 5329.575 rpm c . So that a body force model could be built for an N3-X fan, the R4 fan was scaled up 
to have a fan diameter ratio of 40.25/22 and the blade stagger angle and rotor speed were adjusted to match the 
specifications. The blade stagger and rotor speed were modified, and the Swift code was run to check the modifed 
fan’s performance. Through an iterative procedure, the stagger angle was reduced by 5° and the rotor speed was 
determined to be 5913.667 rpm c . 

In Fig. 19, numerical results from the SWIFT code and the current Euler simulation with the body force model 
are well matched by the body force model. When the body force model are applied here to a viscous flow simualtion 
of the N3-X configuration, entropy production is increased in the viscous simulation because the hub and shroud 
boundary layer effects were duplicated in the Navier-Stokes simulation and the body force model. The body force 
model for Navier-Stokes simulations of N3-X was, therefore, corrected by reducing entropy production by 0.003, 
which was obtained by comparing inviscid and viscous simulations of the R4 fan in the clean inflow conditions. The 
actual amount of entropy correction would vary for different flow conditions and computational meshes. 
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Figure 19. Generation of the body force model of the scaled R4 fan for installation on N3-X. 

B. Flow Simulation Results 

A flow simulation of the N3-X HWB configuration was conducted for a cruise flight condition: Mach 0.84, an 
altitude of 35,000 ft, and an angle of attack of 2°. 

Figure 20 shows sectional side views of total pressure contours for propulsor passages indexed by the order in 
the span wise direction. The region that is swept by the fan blade and has the body force model applied can be easily 
recognized as the steep axial gradient of the total pressure contours around the hub forebody. It can be also noted 
that there is a massive leading edge separation at the nacelle inboard cowl surface and shock-induced flow 
separations in the mid cowl region in the spanwise direction. 

The difference in the ingested boundary layer thickess can be noted from the total pressure contours. The inner 
passages have thicker boundary layers than the outer passages have because of the difference in distances from the 
airframe leading edge. Therefore, the outer passages would have higher fan face recoveries than the inner passages 
have. Figure 21 compares fan face pressure recoveries of the eight passages. As expected, the pressure recovery 
increases in overall as it goes from inboard to outboard. 

Performance of the R4 fan performance was degraded in distorted inflow conditions shown in Section IV B. To 
check the installed performance of the present fan (staggered blade R4 fan) in the mail -slot nacells of N3-X, fan 
pressure ratio and adiabatic efficientcy are presented in Fig. 22. Installed fans behave independently with each other 
depending on the fan face distortion and corrected MFR. 

C. Shape Design Results 

To improve the aerodynamic performance of the mail-slot nacelle, the cowl surface shape was designed to 
minimize the aerodynamic drag using a gradient-based optimization method. Sensitivity analysis is one of the most 
important components to a gradient-based design optimization. The adjoint method allow efficient calculation of 
objective function gradients for a computational cost independent of the number of design variables. In the present 
study, we adopt a discrete unstructured Navier- Strokes adjoint solver 26,27 developed from the flow solver described 
above in Section II. A. The adjoint solver uses GMRES 28 for the time integration with LU-SGS as a preconditioner. 
Parallel processing is also made by domain decomposition of the computational mesh and Message -Passing 
Interface (MPI). 

For geometry parameterization, CST (Class-Shape-Trasformation) method 29 is employed. The CST method 
allows an easy specification of leading edge radius and trailing edge angles of airfoil shapes. The number of design 
parameters was 30: six parameters per a design section times five design sections. 

The overall design procedure is as follows: First, a flow analysis is conducted for the current design 
configuration. Then an adjoint sensitivity analysis is performed based on the flow analysis results to determine a 
search direction. A step size along the search direction is selected by a line search method with the slope along the 
search direction. In the present study, a quadratic polynomial fitting is used for the line search. Then design 
variables are updated using the gradient information and step size. Design surfaces are then modified using the CST 
method, and volume meshes are modified accordingly using a spring analogy. This loop is repeated until the design 
converges, or for a limited number of iterations. As a gradient based optimizer for the unconstrained minimization 
problem, the conjugate gradient method 30 is used. 
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Figure 23 compares Mach contours of initial and design configurations at mid region sections. By the optimal 
shape design, the strengths of shock waves at the mid region on the cowl surface have got remarkably reduced as 
shown in Fig. 23, leading to a drag reduction of 20 counts. But, there still exit flow separations in the inboard region 
as can be seen in Fig. 23 (a). The cowl section shapes became more thicker and convex by the design. More drastical 
shape change of the nacelle such as sweep angle reduction is deemed to be required to remove the large flow 
separation at the inboard region. 




Figure 20. Normalized total pressure contours at center section of each passage 
(Passage index increases from inboard to outboard). 
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Figure 21. Comparison of fan face pressure recoveries 



(b) Fan efficiency vs. MFR 

Figure 22. Performance comparison of clean-inflow and installed fans 
(1: passage 1, 8: passage 8) 
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(a) Passage 3 



(b) Passage 4 
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(c) Passage 5 


Figure 23. Mach contours at center section of selected passages 
for initial(left) and design(right) shapes 
(Passage index increases from inboard to outboard). 
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VII. Concluding Remarks 


High-fidelity flow simulation and optimal shape design were conducted for the N3-X hybrid wing-body 
configuration with the turboelectric distributed propulsion system using 16 fans in the mail-slot propulsor. The mail- 
slot propulsor geometry was generated by an efficient CAD (computer-aided design)-free procedure. A body force 
model generation procedure was suggested for efficient and accurate flow simulations of inlet-fan interaction 
problems. Flow simulation results for the N3-X hybrid wing-body configuration with the body force model show 
the validity of the present approach. Also, an optimal design of the nacelle cowl shape greatly reduced the strength 
of show waves and flow separations on the cowl surface. Inlet shape design for improvement of the fan performance 
will be conducted as a future work. 
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